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Abstract 

We discuss a differential integrable hierarchy, which we call the (iV, M)-th KdV hierarchy, whose 
Lax operator is obtained by properly adding M pseudo-differential terms to the Lax operator of 
the iV-th KdV hierarchy. This new hierarchy contains both the higher KdV hierarchy and multi- 
field representation of KP hierarchy as sub-systems and naturally appears in multi-matrix models. 
The N + 2M — 1 coordinates or fields of this hierarchy satisfy two algebras of compatible Poisson 
brackets which are local and polynomial. Each Poisson structure generate an extended TEi+oo and 
Woo algebra, respectively. We call W(N, M) the generating algebra of the extended algebra. 
This algebra, which corresponds with the second Poisson structure, shares many features of the usual 
Wn algebra. We show that there exist M distinct reductions of the (N, M)-th KdV hierarchy, which 
are obtained by imposing suitable second class constraints. The most drastic reduction corresponds 
to the (N + M)-th KdV hierarchy. Correspondingly the W(N, M) algebra is reduced to the Wn+m 
algebra. We study in detail the dispersionless limit of this hierarchy and the relevant reductions. 
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1 Introduction 



Integrable hierarchies are a topic of increasing importance in theoretical physics. It has been realized 
recently that they play an essential role in the study of 2D quantum gravity and topological field 
theories, as well as in matrix models. In particular matrix models exhibit an extremely rich integrable 
structure. The origin of this interest in integrable hierarchies goes back to 1989 when three different 
groups, through the so-called double scaling limit technique, obtained remarkable non-perturbative 
results in 2D quantum gravity. In particular they found that the partition function of one-hermitean 
matrix model with even potential satisfies the KdV hierarchy equations Later on it has been shown 
that topological field theories coupled to topological gravity and the Kontsevich model also possess 
this integrable structure f2[||- After these successes several authors have conjectured that multi- 
matrix models should be governed by higher KdV hierarchies Q|. Now, it is rather straightforward 
to extract discrete hierarchies from matrix models. The difficulties start when we try to pass to 
differential hierarchies. In particular the double scaling limit technique does not prove as powerful 
and manageable in multi-matrix models as in one-matrix models. 

In H we proposed an alternative approach to investigate one-matrix models, by which we could 
extract a differential integrable hierarchy from the discrete one without reference to any continuum 
limit. Our basic observation is the following one: in one-hermitean matrix model there naturally 
exists a discrete (or lattice) integrable hierarchy - the Toda lattice hierarchy; if we treat the first flow 
parameter as space coordinate, this discrete integrable hierarchy can be re-expressed as a continuum 
(differential) integrable hierarchy, which admits the following Lax pair representation 

L 2 = d + R-^-^, (1.1a) 
^-L 2 = [{U 2 ) + ,L 2 ] (1.1b) 

where d = ^ = -J^-, t\ is the space coordinate, while t r (r > 2) are real 'time' (coupling) parameters. 
The subscript "+" means that we keep only the terms containing non-negative powers of d. R and 
S are the independent fields of the system. This differential integrable hierarchy is referred to as 
two-boson representation of the KP hierarchy^. This hierarchy can be reduced to the KdV hierarchy 
by imposing the second class constraint 5 = 0. 

In ||] we showed that this approach is also applicable to multi-matrix models (with bilinear 
couplings among different matrices) and permits a systematic analysis of these models. Such multi- 
matrix models are characterized by generalized Toda lattice integrable hierarchies which, via the 
same procedure applied to one-matrix models, can be rewritten as the following differential integrable 
hierarchies 

M 1 1 1 

L 2M = d + g *i d _ Sid _ Si _ i ■ ■ ■ (L2a) 

-^-L 2M = [(L r 2M ) + ,L 2M ] (1.2b) 

where a%, om? Si, Sm are independent coordinates (fundamental fields or dynamical variables) 
of the system. These integrable hierarchies are called 2M-field representations of the KP hierarchy. 
It is just these differential integrable hierarchies (together with the string equations) that contain 



*This integrable structure also shows up in WZW model and Conformal Affine Toda field theories (CAT models) Q. For 
a more mathematical approach see |Q and references therein 
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the information of multi-matrix models. In || we showed that they are related to the higher KdV 
hierarchies via Hamiltonian reduction. 

In particular in |J we examined in detail the 4-field representation of the KP hierarchy and 
we obtained two distinct integrable hierarchies by suppressing successively the fields Si , S2 • The 
corresponding Lax operators are 

L21 = d 2 + ai + a 2 n 1 a , L 3 = d 3 + aid + a 2 . 
o - 6 2 

In the present paper we carry our analysis further along the same line. Precisely we will show 
that: 

(i). There exists a general integrable differential hierarchy 

£-L=[(l*).>Q (1-3) 



dt 



r 



with the following Lax operator 

N-l M -. 1 1 

L = d N +^2a l d N - l - 1 + ^2a N+l - 1 - -- — W>1, M>0; (1.4) 

i=1 i=1 o — bio — bi^i a — b\ 

which involves (iV + 2M — 1) independent fields. We call them the coordinates of the hierarchy. Since 



the Lax operator (1.4) is obtained by adding M pseudo-differential terms to the Lax operator of the 
usual N— th KdV hierarchy, and in the lack of a preexisting terminology, we refer to this integrable 
hierarchy as the (N, M)-th KdV hierarchy (or the M-extended iV-th KdV hierarchy) . We refer at 
times to the N + 2M — 1 fields together with all the properties implied by integr ability as 'mode/'. 
So, for example, we will be talking about the (N, M) model in connection with the hierarchy just 
introduced. 

(ii). Secondly we will show that the coordinates of the (N, M)-th KdV hierarchy form two closed 
algebras with respect to the two Poisson brackets which constitute the bi-Hamiltonian structure. They 
generate the extended W\ +00 and Woo algebras, respectively: in other words, there exist combinations 
of the fields and their derivatives which satisfy the extended Wi+oo and Woo algebras. 
(Hi). Finally we will prove that there exist various Hamiltonian reductions. Precisely we show that 
it is possible to suppress the S fields one by one and still obtain integrable hierarchies; at the end of 
this cascade reduction one gets a KdV hierarchy of order N + M. 



The motivation for this work stems from the fact that the pseudo differential operator (1.4) is, after 
reduction to the standard form, the most general operator which appears in multi-matrix models. In 
a companion paper we use the hierarchies of this paper to calculate the correlation functions of 
the corresponding matrix models. 

We remark that there is not only one possible viewpoint to study the system specified by (|1.3|) 
and (|l.4| ). In fact, a posteriori, one will remark that the Lax operator ( |l.4|) can be envisaged as a 



reduction of Lax operator of the type ( 1.2a| ) with M replaced by N + M — 1 and in which N — 1 of 



the S fields have been suppressed. The results of the analysis are of course invariant if we change the 
point of view. 

This paper is organized as follows. In section 2 we present a brief review on the pseudo-differential 
analysis of integrable systems with reference to the general pseudo-differential operator ( |2,1| ) and 



write down the corresponding integrable hierarchy (2.10). In section 3 we will show that this general 



integrable hierarchy ( 2.10 ) admits a particular restriction which leads to the (N, M)-th hierarchy (1.3). 



As stated above, the independent fields of ( |1.4| ) satisfy two (N + 2M — l)-dimensional algebras. We 
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refer in particular to the one corresponding to the second Hamiltonian structure as the W(N,M)- 
algebra. Its properties will be illustrated in section 4. In section 5 we will discuss the reduction 
procedures. We will show that the (N, M)-th KdV hierarchy ( |1.3| ) possesses M different reductions 
which are characterized by M second class constraints Si = 0(1 = 1,2, . . . , M). In particular, when 
we impose the constraint S\ = 0, the (N, M)-th KdV hierarchy and W(N, M) algebra are reduced 
to (N + 1, M — 1)— th KdV hierarchy and (N + 1, M — 1) algebra, respectively. If we suppress all the 
Si fields in succession, we will obtain the following two sequences 

(N, M) - thKdV Sl > =0 ) (N + 1, M - 1) - thKdV S r° > (JV + 2, M - 2) - thKdV 



Sm V 1 ° (N + M -1,1)- thKdV Sa ( ~° > (iV + Af, 0) = (N + M) - thKdV. (1.5) 
for the hierarchies, and 

W(N, M) S T° > W(iV + 1, M - 1) S2 > =0 > W(JV + 2, M — 2) 



^"^-> W(AT + M - 1, 1) W{N + M, 0) = W N+M . (1.6) 

for the W(iV, iW) algebras, respectively. The double arrow — ► — * means reduces to throughout the 
paper. 

Our results on reductions can be nicely summarized by means of a Drinfeld-Sokholov representa- 
tion, section 6. 

We also consider the dispersionless version of the (N, M)-th KdV hierarchy in section 7 and 8. 
They constitutes the genus part of the hierarchy in the context of matrix models, sec |l(|. We 
find that the dispersionless (N, M)— th KdV hierarchy (|7.2j ) admits, in addition to the ones already 
considered, another type of reductions, in which some Si fields are identified with one another instead 
of being suppressed. In this way we get more possible reductions. All of them can be seen as 
restrictions, i.e. we can obtain the reduced hierarchies by imposing the constraints on the Lax pair. 
Section 9 contains some final remarks. 



2 The pseudo— differential analysis 

In this section we give a short review of integrable differential hierarchies, in particular we collect 
some results we will need later on concerning pseudo-differential analysis, coadjoint orbit method and 
bi-Hamiltonian structure. We refer to existing books and reviews (see [|ll]] and references therein) for 
more detailed explanations. 

2.1 Pseudo-differential analysis 

Let us begin with the most general pseudo-differential operator (denoted PDO) 

N 

A = Y J <x)8 i . (2.1) 
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where Ui(x)'s are ordinary functions. As usual d = J| is the derivative with respect to the space 
coordinate x, while d~ l is a formal integration operator, which has the following properties 

d~ 1 d = dd- 1 = 1, 

d- j - 1 u = jr{-iy r r + w juWa-j'-"- 1 , (2.2) 



where u' = || , u" 



<9 2 « 

9a? ' 



i)=0 \ / 

, = The above formula together with the usual Leibnitz rule, 

du(x) = u(x)d + u'(x), [d, x] = 1 



provides an algebraic structure on the whole operator space formed by the general PDO's (|2.1|). We 
call this algebra pseudo-differential algebra p. 

For any pseudo-differential operator A of type ( |2.1| ), we call order the number N (the highest 
power of d). We will use the conventions 

A + : pure differential part of A; 
A_ : pure integration part of A; 



.4 



« > i; 



.4 



(0 



-oo < i < N. 



Therefore any pseudo-differential operator has a natural decomposition 

A = A + + A_. 
which leads to an analogous decomposition in p 

P= P+ + P- 



In terms of this truncation, we can introduce the useful mapping []15| 

TLA = A+-A-. 
This mapping defines another algebraic structure on p 

[X,Y] n = l -([KX,Y] + [X,KY])= [X + ,Y+] - [X_,YL], VX,y G p. 



(2.3) 



(2.4) 



This 7?.-commutator will be very important in our discussion on Hamiltonian structures. But before 
that, we need more notations. We call n_i(x) the residue of pseudo-differential operator A of type 
(|2.1|) and denote it by 



res^A = u_i(x) or 
The following functional integral 



A 



(-i)' 



< A >= Tr(A) = J (resA) = / u_i(x)dx 



(2.5) 
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naturally defines an inner product on the algebra p, which is nondegenerate, symmetric and invariant 
Using this inner product, we may express a functional of the fields Uj's as follows 

f x (A)=<AX>=A(X), X = J2d i Xi (x) 

i 

where Xi( x )' s are testing functions, and X may be thought of as a cotangent vector at the point A. 
Generally, the cotangent vector df is defined as 

5f(A) = f{A + 5 A) - f(A) d ^< df, 5A > . 
We denote by J-{p) the functional space formed by all functionals fx defined as above. 

2.2 Bi— Hamiltonian structure 

By means of the 7£-commutator introduced in the previous subsection, we may define the adjoint 
action of the algebra p on itself^ 

AdyX^ [X,Y] n . 

The coadjoint action of p on the functional space .F(p) is specified by 

Ad* Y f x (A) d ^ f A(Ad. Y X) = A([X,Y] n ), 

For a fixed fx, as Y varies in p, Adyfx spans an orbit in the functional space J-(p) which is called 
the co-adjoint orbit. Since we may view Y as a co-tangent vector, this co-adjoint action naturally 
defines a Poisson structure on T(p) 

{fx,fYh(A) = A([X,Y] n ). (2.6) 

With respect to this Poisson bracket, the conserved quantities (Hamiltonians) have very simple form 

N 

H r = — < An >, Vr > 1; (2.7) 
r 

The corresponding cotangent vector at the point A reads 

r-N 

dH r = A (1 "_ N>oo) , Vr>l. (2.8) 
The time evolution of a function fx (A) is given by 

^fx(A) = {/x(A),# r+JV }i =< A, [X,dH r+N ] n >, 



*The product is symmetric since < AB >=< BA >, while invariance means < A[B,C] >—< [A,B]C >. Occasionally 
we will also denote the product in other ways: < AB >= A(B) — Tr(AB). 

t The usual adjoint action is AdyX = [X, Y], here we use the same notation to denote the adjoint action generated by 
the 7?.-commutator. 
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From eqs.(|2.4 2JS), we see that 

< A, [X,dH r+N ] n >=< A, [X+, (A*)+] - [X^A*^^] > 

= -< [A,(A7r) + ],X + > + < [AAf^^X- > 

= -<[A,(A%) + ],X + > + < [A,(A*)-.],X- > 
=<X,[{Att)+,A] > . 

In the second and the last equalities we have used the invariance of the inner product; in the third 

r 

step we have added a vanishing term < [A, AW^ _j~],X- >, since 

[A^ ( * 00 ,_ iV) ]ep-, and <Y,Z>=0, if Y,Zep_. 

Therefore we have 

^f x (A)=<X,[(A*) + ,A}>, (2.9) 

Suppose that X is time independent; then we obtain the time evolution equations of the pseudo- 
differential operator ( |2.l[) 

£-A=[(A*)+,A], Vr>l. (2.10) 

This is the general integrable differential hierarchy we have obtained from pseudo-differential analysis^. 
In order to prove its integrability we derive the second Hamiltonian structure, |l3|], compatible with 
the first, i.e. a second Poisson bracket such that 

{/X,#r+Ar}l = {/X,#r}2. (2.11) 

The LHS is just (|2.9|). Taking the residue of the following equality 



we get 



An) ,A] = [A,[An i 

/(l-JV.oo) 1 1 \ J (-00 -N) 1 



N(A i - N) ) = [A-_ N ^ ) ,A] { _ 1) . 



On the other hand, 



r r — N / r~N . , — jy , 

/ r — N \ / r—N\ 

[(A—) {1 „ N>oo) A)+(A—) 



, AT s 

'(-AT) 



*If we define another kind of "/^.-operator such as 

IZlX — Xi <oa — X-aaj-i, 

we obtain non-standard integrable hierarchies for some values of Z@. 
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Therefore, we have 



LHS of eq.(2.11) =< AX(A*)+ > - < XA(An) + > 



=< AX 



f / r-N 

I A[A— 



(1-Af,oo) / + 



> - < XA A 



> + <(ax-xa)(a^) ( _ n > 



=< Ax(AdH r ) + > - < XA(dH r A) + > +— J (VV#r-,^](-i) 

=< (XA)+dH r A >-< {AX) + AdH r > +1 / [A dtf r ] ( _ 1} (^[A *](-!) 

Now we can see that all the terms in the last equality are linear in the cotangent vector dH r and 
in X. This is therefore a candidate for a second Poisson bracket. If we replace dH r by an arbitrary 
cotangent vector Y, we have the general expression of it ( 2.1C| ) 



{fx,fyh(A) 



= < {XA) + YA > - < (AX) + AY > 

+ i J^YU^d-^xu 



(2.12) 



One can show that this Poisson bracket satisfies the Jacobi identity [14], so it is indeed a well-defined 
Poisson structure. With respect to this bracket, the conserved quantities are the same as in eq.(p.7|), 
and generate all the flows in the hierarchy ( |2.10 ). Therefore we proved that the system ( |2.10| ) possesses 
a bi-Hamiltonian structure and is therefore integrable. 



3 The (TV, M)-th KdV hierarchy 

In the previous section we have constructed the bi-Hamiltonian structure for a general pseudo- 
differential operator and the corresponding integrable hierarchy. In this section we will prove that 
the general integrable hierarchy ( |2.10| ) admits a particular restriction which leads to the (N, M)-th 
KdV hierarchy Ql.3| ) with the Lax operator (|1.4|). 



3.1 The consistency of the general flows 



In order to show that the restricted Lax operator ( |1.4| ) preserves the hierarchy ( 2.10| ), and gives in 



fact the (N, M)-th KdV hierarchy, we should first prove the consistency of the flows defined by ( |1.3| ), 
and then check that they commute. 

For the usual iV-th KdV hierarchy ( KP hierarchy), the scalar Lax operator is a pure differential 
operator ( pseudo-differential operator), so it is very easy to see its form invariance under the time 
evolutions. For instance, the scalar Lax operator involved in the iV-th KdV hierarchy is 

N-l 

N-l-1 



A NK dv = d N aid 
i=i 



Obviously we have 



r 

[(^NKdv) N ; ^NKdv] = [^-NKdV i (^NKdv) 
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The LHS is a purely differential operator, while the RHS indicates that this operator is of order N — 2. 
This enables us to define the flows as follows 



TTT-^NKdV = [(^NKdv)^>^-NKdv]- 

Therefore, the consistency of the flows in the iV-th KdV hierarchy (or KP hierarchy ) case is very 
simple. 

However the present case is much more complicated since the scalar Lax operator ( |1.4j) contains 
both differential and pseudo-differential parts. It is a nontrivial result that the time evolutions of the 
Lax operator (|1.4|) are form invariant. This subsection is devoted to proving exactly this. 



Proposition 3.1 The LHS and the RHS of eq.(1.3) are compatible. 

Proof. To start with we write down the r-th time evolution of the Lax operator ( |1.4| ) 

d_ y^/dai\ b n-i-i f do-N+i-i \ 1 1 1 

dt r fr[^ dt r^ i^l dtr ' 9-Sid-Si-i'" d-St 

+ l^^- 1 ^---^©^---^ (3 - 1} 

The problem is to prove that the LHS of eq. (|1.3j ) contains the same type of terms. Our first task will 
thus be to classify the terms in RHS of this equation. To this end we introduce a new basis set for 
the pseudo-differential operator algebra p. We split the basis set into several classes: 
the differential operator class: we adopt the usual basis 

{f i d i ,ieZ+}, 

with arbitrary functions /j's. By definition this basis spans the operator space F°. 
the first class of integration operators: 

where f{s are arbitrary functions. For a fixed I, the above operators span by definition the space i*} . 
We also define the direct sum space F 1 = 0/i 7 ) 1 . 
the second class of integration operators: 

flk-^r- — . . . — — - — Ot— ^ — . . . — — - — , KKM: 1< k < I. 
Jlk d-S l d-S k y d-S k d-S^ - - ' - - 

Once again /^'s are arbitrary functions. We denote by Ff k the space spanned by the above operators 
with fixed l,k, and 

F 2 = {®i, k F£ k , 1 < I < M; 1 < k < I}. 

All these operators are linearly independent. Actually they do not constitute a complete basis of the 
algebra p. However, since the residual basis subset will not show up in our later discussion, we do not 
need to write them out explicitly. Going back to eq.( |3.1[ ), we immediately see that its RHS contains 
all these terms. We will say that a pseudodifferential operator takes the standard form F°, F 1 or F 2 
if it belongs to the corresponding vector subspace. 
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What we should do next is to prove that the RHS of eq.([ll|) can be recast into these standard 
forms. 

In order to simplify the calculation of the commutator in eq. ([l]^) , we adopt the notation 



r-2 



0{r) H f (L^) =^ + 5>^, Vr>l. 



1=1 



where the a/'s are certain functions of the {a;}'s and the {S/}'s. The commutator in eq.(1.3) can be 
decomposed into two parts 

[0 {r) ,L] = [0 {r) ,L + ] + [0 {r) ,L„}, 

Obviously the first part is a purely differential operator of order N — 1, while the second part can be 
further simplified 



M 



[0( r ),L_] = [Q( r ) , gjy 



1 



1 



1 



l=i 
1 



d-Sid-Si-i "' d-Si 
1 1 



+Z-1J 



i=i 

M I 

EE 

1=1 k=l 



d-Sid-Si-x d-S 



aN+i-i 



d - Si ■ • • d - s k 



r o (r) ,d-s k ] 1 



d-Sk'^d-Sx' 



(3.2) 



Here we have used 



1 -[o [r) ,d-s k ]- 1 



d-Sk 1 d-s k l ^' ^d-s k 

The remaining two commutators in the expression (|3.2|) only involve purely differential operators, 



both of them are of order r — 1. Therefore we find that the RHS of eq.(L3) contains the following 
three kinds of terms 



r g. 



o • 



purely differential operator part [0( r ),L H 



j(r) . r-) J i. i 



(r _i 



(0) d-Si ■ ■ ■ d-S k 0(r—l) d-S k ■ ■ ■ d-S! ' 



1 < I < M: 



1< k < I < M. 



where 0( r ) is an arbitrary purely differential operator of r-th order, while O(o) is an ordinary function. 
We will denote 

Gi = {]>>, 1</<M}, 
i 

G , 2 = {E J / { !T 1) ' i<k<i<M}. 

i,k 

We will need one more operator space, G 3 , which contains the following type of terms 



V 



(r-l). n 1 1 n 1 1 



l.k 



1 < k < I < M. 
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Here Or r -\) is another purely differential operator. In the case k = I, this is the same as Gi, i.e. 

Vff=l! r) . (3.3) 

In the remaining part of this subsection we will show how to recast the above terms into the standard 
forms F°,F\F 2 . 

Go terms: 

We notice that 

[(i*) + ,L] = [L, (!*)_], 

so the power expansion of the above commutator has order < N — 1, i.e. Go-terms are pure differential 
operators of (N — 2)-th order, i.e. they have exactly the same form as .F°-terms. Therefore Go-terms 
contribute to the time evolutions of the fields ai (1 < I < N). 

G\ terms: 

In order to simplify Gi-terms, we recall that for any purely differential operator of r-th order, 
and any ordinary function /, we can find two (r — l)-th order purely differential operators C( r _i) and 
Ou.-i) such that 

{r) = {d-f)0 {r _ x) +g, and (r) = 6 {r _ x) {d - f) + g. (3.4) 

where g and g are two ordinary functions specified by the above equalities. The second equality 
immediately leads to the following results 

repeating this procedure we can reduce all the Gi-terms to F 1 (possibly F°) terms. For instance 

[ + , a N+l ^} = (r _ 1} =/ + O (r _ 2) (d-Si) 
= h + (fx + (r . 3) (d - 5,-0) (8 - (3.5) 



l-l 

= /o + E M 9 ~ Si-i+i)(d - Si-i) ...(d-St) 
i=i 

+0 {M) (d-S 1 )(d-S 2 )...(d-S l ), 

where all the functions fi's and the operator Or r _i_i\ are completely determined by the commutator. 
Therefore the Gi-term 

UT^ 1 1 1 1 



'd-Sid-Si-i d-s x 



becomes an F — term. If r > I + 1, the last term in eq.( |3.5| ) contributes an F -term Or r -i—\)- We 
may summarize the procedure with the following diagram 

T (r-I-1) T (r-0 

(3.6) 



r(r) 
L l 


A Z-1 


T ( r " 2 ) ► 


T (r-I-1) 


T (r-0 


I 


1 


I 


I 






FU 


^-2 
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Go terms: 



Now let us consider the G^-terms. Their form is 



,(r-l) 1 
J Lk = a N+l-l 



d-Si d-s k 



' [(L$)+,d-S k ] ' 



d-S k '"d-S 1 



OjV+Z-1 



a-N+i-i 



d-sr" d-s k 0( - r - 1) a-s k "'d-s 1 



i 



-o 



1 



a- Si'" a- s k ~ {oy a -s k '" a-s 



+o>N+i-r 



1 



-O 



1 



1 



■a-Si'" a-Sk+r^a-Sk'" a-Sx' 

in the last step we have used the first equality in eqs.( |3.4j ). This decomposition simply means 



r(r-l) 

'Lk 



V 



(r-2) 
Lk 



(3.7) 



Here and in the following ==> means decomposition. In other words, G^-terms can be decomposed 
into F 2 and G3 terms. Therefore it remains for us to treat G^-terms. 

Ga— terms: 



As we know, a Gr3-term is of the form 



1 



-o, 



1 



i,k - "N+i^dZTs; ■ ■ ■ a - s k+1 ^ r -^a^s~ k ■ ■ ■ s^sT 

With a simple calculation we get 



1 



a — s k +i 



C(r-2) — C(r-2) 



1 



a-s k+1 L d-s, 



k+l 



-,o (r _ 2) ] 



o 



1 



(r-2) 



a — s k +i a — s k+ 



[d - S k+ x, 0( r -2)] 



o 



1 



1 



1 



(r-2) 



a — s k +i a — s k+ i a — s k+ 



■Ou. 



a — s k+ i 
1 



(r-3) 



a — s k+ i 



This result shows 



V 



(r-l) 
Lk 



v (r-2) T (r-1) 
v l,k+l w J Z,fc+l ' 



(3i 



The crucial feature of this step is that we have moved the operators in G3 one step to the right. 

Combining the procedures (|3.7|) and (|3.8|) , and remembering the fact (|3.3[) , we can recast the 
G2, G^-terms into the standard forms. Diagrammatically 



V (r-1) _ 


v (r ~ 2) - 

^ v I, k+l ~ 


v (r ~ 3) 


( r _Z+fe_l) 


T 


\ T 


\ T \ ■ 


■ \ T 


j(r) 
J l,k 


T (r-1) 
J Z,/«+l 


T (r-2) 
J Z,fc+2 


T (r-Z+fc) 


i 


I 


1 


1 


F l,k 




P 2 

r /,fc+2 


F 5 



r (r-Z+/c-l) 
l Z 



(3.9) 
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The last term in the first line can be treated through the procedure (|3.6| ), So we finally proved that 
the RHS of eq.( |1.3| ) has exactly the same form as its LHS. Comparing their explicit expressions, we 
can obtain the equations of motion of the fundamental fields. 



3.2 Commutativity of the flows 

Our next problem is to prove the following 

Proposition 3.2 All the flows defined in (11.31) commute with one another, i.e. 



Proof. This is elementary to prove since we have explicitly 



RHS = [[L$ ,Ln] + ,L] + [L» ,[Lf ,L]\ 

r I r _Z_ _l_ 

4])-[L,[L%,„ + 



[[LH,L"} + ,L] - [L»,[L,L»]] - [L, [L»,L*" 



= [[L$,LV] + ,L] + [L$,[L$,L]] = LHS. 
In the second equality we have used the Jacobi identity. 

3.3 The non— linear evolution equations 

Now let us give here the first non-trivial flow of the hierarchy ( |3.11a| ) (see Appendix A for the 
derivation) , 



d n 2 N 

— a l = ai +2a l+1 --l l + 1 



I e ( N r_ k 1 ) afc °^ fc) ' 1 ^ ^ n - i; (3 - na) 



fc=i 







l+l 



—a N+ l = a N+l + 2a' N+l+1 + 2a' N+l Si+i + 2a N+ i(Y\ ' < / < M — 1; (3.11b) 
dt2 fc=i J 



3 2 

— S l = -a' 1 + 2S l S' l -S'i -2(Y,S k ) , 1 < I < M. (3.11c) 
dt 2 N V fc=1 / 

This set of non-linear evolution equations are crucial objects, all the important ingredients like the 
bi-Hamiltonian structure and the algebraic structures, as well as the conserved quantities, are rooted 
in these equations; the pseudo-differential analysis is only a powerful tool to make them explicit. In 
other words, we can say that the above equations define a bi-Hamiltonian system with two Poisson 
structures (p^) and fl2.12|) , the hierarchical equations (|1.3|) can be thought of as symmetries of these 
non-linear equations. 

We end this section with a few remarks. 
(i). Comparing the first — 1 equations ( [3.11a| ) with the A-th KdV equations, we see that they are 



i-i 



the same, except that the on field is involved in the time evolution of the field a^-i in eq.( 3.11a ). 



Therefore, if we set = Si = 0,1 = 1,2, ... , M; we recover the A-th KdV hierarchy. So we may 
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view this set of equations ( 3,lla|) — ( 3.1 lcj ) as a generalization or a perturbation of the iV-th KdV 
equations by means of the fields {ajv+z-i, Si] 1 < I < M}). 

(ii). In the N = 1 case, the integrable hierarchy ( |1.3| ) is nothing but the 2M-field representation of 
KP hierarchy || . 

In this sense we can say that the (N, M)-th KdV hierarchy contains both A^-th KdV 

hierarchy and the multi-boson representations of KP hierarchy. 



4 W(N, M)— algebra (or the extended W— algebra) 

As we know, the second Hamiltonian structure of the A^-th KdV hierarchy gives rise to a Wn algebra, 
and the second Hamiltonian structure of the KP hierarchy leads to a Wqq algebra. We will see that the 
bi-Hamiltonian structure of the integrable hierarchy ( |1.3| ) results in two finite dimensional algebras, 
which generate the extended W\ +(X3 and algebras, respectively. 



Before we proceed we wish to clarify the meaning of our choice of coordinates in (1.4). Up to 
now in fact we have only used {a/, I = 1,2, . . . , N — 1; cijv+2-i, Si, I = 1,2, ... , M} as our dynamical 
fields (coordinates). However, one might choose any other system of coordinates. For instance, we 
can introduce a new set of coordinates in the following way 



Noting 



(Inn)' = -Si, qi = ^^, 1<1<M; (4.1) 

n 



l l 

;n = n- 



d d + (In r^y' 

we can immediately rewrite (Fj) as follows 



N-l M 



L = d N +J2 aid N ~ l ~ l + J2 Qid- 1 (— ) d- 1 . . . p) d- l n . (4.2) 
i=i i=i ri ~ 1 ri 

Each set of coordinates have their own advantages (and in general different physical meaning). For 
example, the first choice (|1.4| ) leads to simple and local Poisson algebras, while the Poisson algebras 
in the second choice ( |4.2[ ) will contain non-local terms. We can think of the passage from one set 
of coordinates to another set as a field-dependent gauge transformation of the corresponding linear 
system (see the representations of the Lax operators in terms of linear systems of section 6). For this 
reason we will refer to different choices of the coordinates as "different gauges" . 
We notice however that we can write 

JV-l oo 

L d ^ 8 N +Y, a^- 1 - 1 + £ u n d- n -\ (4.3) 



where 



1 = 1 71=0 



^ = E E a N+ i^(-d+sir(-d+si^r-- ...(-d+s^ -i, 

1=1 £!=iw=»-i+i 



< n < M - 1: 
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M 



II 



n = E E a N+l - 1 (-d + S l r>(-d + S l „ 1 r^...(-d + S 1 )^-l, (4.4) 
1=1 ELi w^-'+i 



M < n < oo. 



The /orm of the operator L is gauge invariant and by (a suggestive) abuse of language we call the 
coordinates {a;, 1 < I < N — 1; ui, I > 0} 'gauge invariant'. 

The expressions ( |4.4| ) have the following remarkable features: 

• the wj's are linear in the a/ (1 < / < N + M — 1) fields, but highly non-linear in the Si fields; 

• the iti's do not contain any derivative of the a/ (1 < I < N + M — 1) fields while they contain 
derivatives of the Si fields. 

• these formulas show that the subsets of fields {a/, 1 < I < N — 1} and {oat+;, < I < M — 1} 



introduced in (1.4) are of quite different nature, since the former set is canonical the other is 
not. However for later convenience we will keep the notation {ai} for both subsets. 

• in the gauge ( [4.2|) , the u;'s will be linear in the qi fields, but non-linear and non-polynomial in 
the ri fields. 

In the following we will mostly work with the gauge ( |l.4[) . 



4.1 The extended algebras 



If we substitute the Lax operator (O) into (|2.6| ) and ( 2.12 ), we will get two infinite dimensional 
algebras, which we will call extended Woo algebras. The calculation is in principle straightforward, 
but it requires some skillful use of suitable techniques to simplify the results. However we will omit 
any detail and only give the results. 

The extended Wi +oc algebra: 

Proposition 4-1 The first Poisson structure (|2.6|) leads to the following explicit algebra 



{ai,cij}i 



{a i ,u j } 1 

{Ui,Uj}l 



N-l 



y+j-N 

i-l-1 



1 + f— IV 



Qi+j-N+l 



7V-1 



+ E 



1=3 



j-l-l 



a, 



Qi+j-l-N 



5(x - y) 



[± 

i=i 



SVh + EC- 1 )'" 1 / )u l+3 -id l ]5(x-y). 



i=i 



(4.5a) 



(4.5b) 
(4.5c) 



In the above Poisson brackets (as well as in subsequent ones) we use a shorthand notation. They 
must be understood in the following way 



{f(x),g(y)} 1 =0(x)5(x-y), V/, 5 ; 
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Moreover either the indices of ai (ui) fields are in the region (1 < I < N — 1) (7 > 0) or the 
corresponding terms are understood to be absent. Similarly, either the powers of d are non-negative 
or the corresponding terms are absent. We remark that a central term is contained only in eq.( 4.5a ) 
The above Poisson algebra is a direct sum of two sub-algebras, the sub-algebra Q4.5a ) coincides 
exactly with the first Poisson algebra in the 7V-th KdV hierarchy. From eq.( 4.5c ), we see that 

{ui, «i}i = (ui«9 + du 1 )5{x - y), 

Therefore u\ can be viewed as conformal tensor of weight 2. The Poisson bracket {u\,uq}i indicates 
that uq has conformal spin 1. Such sub-algebra ([Ok]) is commonly referred to as a W^i+oo algebra. 

The extended algebra: 



The second Poisson algebra is much more complicated than the former one. 
Proposition The explicit form of the second Poisson brackets is 



{ai,aj} 2 



<•,'>'•'• 1 • E (-i: 



i+i 



l=i+j-N+l 



+ E h ljl) a i+3-l& 



i+j-N 

+ E 

i=i 



d l u, 



+j-l-N + ( — ( l I u i+j-l-N 



i+j-1 



+ 



/ 1 J ~ 



i+l+1 



l=i+l 



N \ I N + l-i-j-1 
i + 1 j I l-i 



i+j-i 

+ Yl h li) dla i+3-i 

l=i+j-N+l 



(4.6a) 



i—l , i — l — 1 

+E(E(-i: 
i=i v k=i 



k+l 



i-l-1 
k 



i+j-l-2 

a i+j~i-k~id h ai + ^2 tf,jik a id k a i+j-l-k- 
k=i+j-l-N 



i-l i+j-l-N-1 

+E E ((-i) fc+1 

1=1 k=l 



k 



u i+j -i^ k -N-id k ai + 



j-l-1 
k 



1 



j+k+ 



! / N-l-1 \ / iV-Jfe-1 



1=1 k=l 



i - I 



j - k 



a, 



Qi+j-l-k-l 



a k 



aid u i+ j^i^k-N-i 
S(x - y); 



{ai,Uj} 2 



N+j 

E 

i=i 



N t 3 )d i ui+3-i+Yi-i) i+i [i j^i-^ 



N-l N+j-l-1 

E E 

i=i k=i 



N + j -l-l 
k 



1 



tyd Ui+j-l-k-l - T7 E 
i=0 



i + 1 M J 



J \d i+j ~ l ui 



i-l i-l-1 



+E E (-D 



fc+i 



i=i fc=i 

i i-lj-l 

--YY 



i-l-1 

k 



Ui+j-i-k-id k ai 



1=1 k=0 



N-l-1 
i-l 



J 
k 



a, 



Qi+j-l-k-l 



Uk 



S(x - y); 



(4.6b) 
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{Ui,Uj} 2 



N+j 

E 

1=1 



N+j W« 

I C U i+ j +N _i 



+ £(-W - v + / 



+E E 

i=0 V fc=l V 
AT-1 JV+j— I— 1 

+ E E 

i=l fc=l 



i - 1 - 1 



E^ 

1=1 

i-l-1 



I I u i+j+N-l 



0' 



uid k u i+j ^- k ^i + £ (-1) 



k+l 



k=l 



N + j -l-l 



N-l N+i-l-1 

+ E E (-i) fc+1 
i=i fe=i 



/ iV + i - z - 1 v 



i=i fc=i 



I \ k 



Ui+j+N-i-k-id k ai 
5(x - y). 



j-k 



where 



Cij N \ i + 



Ui+j+N-l-k-2d Ul 



{ - 1)>( \)+±(-^{ N l ){ N N + r j ' 



(4.6c) 



,1 = ( N \ N + l-i-j-l 

ijl N j + 1 Z-j 



H»-i 

+ E (-i) 

k=0 



k+l N + l-i-j-l \ N + l-j-k-l 
\ k \ N-j-1 



E(-!) 



k=0 



l+k N) N + l-j-k-l 
[ k \ N-j-l 



i-l-1 

tfjik = (~i) 

r=0 



r+fc / iV - Z - 1 \ / iV + Jfe-i-r-1 
I r [ N-j-l 



From eg, (4.6a) can immediately extract a Virasoro subalgebra 



{a i . in } 2 - ( - ( 1 ] d 3 + aid + dai )<)(.r - //). 



(4.7) 



i.e. ai can be interpreted as a semi-classical energy momentum tensor. The Poisson brackets between 
ai and the other fields tell us the conformal dimensions (or spin contents) [ • ] of our coordinates 

Mconf = I + 1, i = l,2,...,JV + M-l; 

N conf = l, Z = 1,2,...,M; 

Mctmf = N + l + l, I = 0, 1, . . . , oo. 

If we assign to ^-function a conformal weight 1, then we have 

[{) }2]conf = 0. 
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The algebra ( 4,6a - Ok] ) contains fields with spins from 2 to infinity, and is linear or bilinear in the 
gauge invariant functions. There exist central extensions represented by the coefficients Cij. For this 
reason we call this algebra the extended algebra. 



4.2 The finite dimensional algebras associated to the (N, 1)— th hier- 
archy 

The above algebras are independent of the particular coordinatization we choose for the Lax operator. 
However, the physical meaning of a hierarchy may essentially depend on the gauge, in particular on 
the number of fields. Therefore we are very much interested in the algebras formed by the coordinates 
we choose. In particular the independent fundamental coordinate fields in the integrable hierarchy 
(|1.3| ) are finite. Therefore we expect two finite algebras corresponding to the two compatible Poisson 
structures. They in turn generate the infinite dimensional algebras of the previous subsection. We 
can derive these algebras from the infinite dimensional algebras ( |4.5a| - 4.5c ) and (|4.6a - 4.6c ) by making 
use of the expressions ( [4.41) . As a example, we consider the (JV, l)-th hierarchy, in which the scalar 
Lax operator is 



N-l 



L = d N +Y.aid N - [ - l +a N 



1 



l=i 



d-Sx 



(4i 



The gauge invariant functions are {a±, a 2 , ■ ■ ■ , a^v-i} and a set of infinite many functions( they are 
generated by only two fields ) 



ui = a N ai, ai = (-d + Si) 1 ■ 1, l>0. 

Proposition 4-3 The first Hamiltonian structure leads to the following Poisson algebra 

{ai,aj}i = (4.5a), l<i,i<JV-l; 
{ai,a N }i = 0, i = 1,2, . . . ,N; 
{ai, Si}i = 5i, N 5'(x -y), i = 1, 2, . . . , N; 
{Si, 5i}i = 0, 



(4.9) 



(4.10a) 
(4.10b) 
(4.10c) 
(4.10d) 



This (iV+1) dimensional algebra generates the W^i+oo algebra ( 4.5a - 4.5c ) through the transformation 
Proposition 4-4 The second Hamiltonian structure leads to the following Poisson algebra 



{ai,aj} 2 = (4.6a), i,j = l,2,...,N— 1; 
{a N ,a N } 2 = {u ,u }2; 
{ai,a N } 2 = {ai,u } 2 , 1 < % < N - 1; 
N + 1 

{Si,Si} 2 = ^—5'(x-y). 



{ a j,Si}^ 



j'-i 

1=0 



-I 



N 

j_( N + l 
N\ j + 1 

N 

3-1-1 



(4.11a) 
(4.11b) 
(4.11c) 

(4.11d) 



d j+i + Y^(_ N + 1 )j- N ( l + 1 ) ( N-l 
i=i 



N(N - I) 



N-j 



aid 3 ' 



aiSi 



N + l 
J -I 



Q'i 



O+B-i: 

1=0 



aiS[ j ~ l) (4.11e) 
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Z=l fc=0 



*(z-y), i<i<JV-l; 



JV-l 

{a*,Si} 2 = + + £ ajCS + ^^-'-^'Ca-y). (4.11f) 

i=i 

The proofs of the above two propositions are not very difficult, so we skip them. 

4.3 The general W(N, M)-algebra 

The Poisson algebras defined by the first Hamiltonian structure is relatively simple. On the other 
hand the second Hamiltonian structure certainly plays a more important role (see the Hamiltonian 
reduction below). Hereafter we only pay attention to the second Hamiltonian structure, and as 
we mentioned before, we denote by W(N, M) the finite algebra represented by the second Poisson 
brackets of the fundamental fields of the (N,M) model which is encoded in (4.6a- Pk| ), Contrary 



to the latter the W(N, M) algebras are gauge-dependent, i.e. they depend on the coordinates we 
choose for the model. In Appendix B we give several explicit examples of the simplest W(N, M) 
algebras. In principle we can calculate any W(N, M) algebra. Unfortunately we cannot exhibit a 
compact explicit form of W(N, M) with arbitrary N and M. However it is not difficult to extract 
some general properties of these algebras. This is the aim of the present subsection. 

{%). W(N,M) in the gauge ( |Q]) : 

In this gauge 

• W(N, M) algebras are local and polynomial , i.e. the Poisson brackets contain neither integration 
operators <9 _1 nor fractional or negative powers of the coordinates. 

• W(N, M)-algebras are linear or bilinear in the a; fields but, in general, highly non-linear in the 
Si fields. 

• The fields are characterized by their conformal spin; there are (N + M — 1) fields with spin 
ranging from 2 to (N + M); in addition, there are M spin one fields. 

• For any Poisson bracket, the <9°-terms in its RHS are either absent or contain derivatives of the 
fields. In other words, for any two coordinates / and g, if 

{/,<7} 2 =B5(x-y), B = J2kd l 

then either bo = or bo contains derivatives of the fundamental coordinates. 

The first and the third assertions can be checked case by case, and in fact they are true for all our 
examples in Appendix B and the previous subsection. The second property can be obtained from the 
transformation (4.4) and the algebra (4.6a-4.6c). Since the algebra ( |4.6a - 4~6c| ) is at most bilinear in 



the gauge invariant functions, while the ui fields are linear in the a\ fields and non-linear in the S\ 
fields, then W(N, M) must be at most bilinear in the ai fields but, in general highly non-linear in the 
Si fields. 
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(ii). W{N,M) in the gauge (^2|): 

• W(iV, M) algebras are, in general, non-local and non-polynomial. 

• Spin content: there are (iV — 1) fields with integer spin running from 2 to N, and 2M fields with 



spins taking value 



\n 1 — "+ 1 

WJconf — 2 ' \ < 1 < M 

hlconf = / + 1 < Z < M. 



2 ' 

In Appendix C we will give some explicit examples of algebras in this gauge, 
(ra). W(N, M) algebras and Wn~ algebras 

W(N, M)-algebras are intimately related to Wn— algebras in the following sense: 

• They share the same Virasoro sub-algebra 

• The fields in the W)v — algebras are gauge invariant, while the W{N, M) algebras contain all these 
gauge invariant fields plus some gauge dependent fields, whose spins are integers or half-integers 
(depending on the gauge choice). 

• (|4.6a| ) is exactly a Wjv-algebra modulo u^-dependent terms. 

• In the next section we will show that the W(N, M) algebra can be reduced to the usual Wn+m 
algebra. 

(iv) . Recurrences among W(N,M) algebras 

From Appendix B, we find that W(l, 2) is almost a sub-algebra of W(l, 3), the only discrepancy 
being that {02,02)2 in the latter case contains a3-linearly-dependent terms. This property holds for 
two more pairs: W(2, 2) and W(2, 1), W(3, 1) and W3. In fact, in general, we can write 

W N cW{N, 1)C . . . CW{N, M)cW(N, M). 

This is to be interpreted in the following way: for any M < M, let us pick out from W(N, M) the 
Poisson brackets among (a\, . . . , a N , Si, S2, ■ ■ ■ , Sj^), and mode out (a N , ^, . . . , ajy+M-i! 
*^Af+2' • • • ' 5 , M) _ dependent terms; then we get the W(N, M)-algebra. 

(v) . M dimensional subalgebras 

Beside the Virasoro algebra fl4.7|), W(N, M) has another very simple subalgebra, the Poisson 
algebra formed by the Si fields 

{Si,S j } 2 = (6 ij + jj)6'(x-y). (4.12) 

To derive it we notice that on a simple basis of dimensional counting the only possible brackets are 

{Si, Sj}2 = const. 5'(x — y), (4-13) 

Then such form must remain unchanged when we take the dispersionless limit (see below). In section 
7, when we analyse the dispersionless {N, M)-th KdV hierarchy, we will prove that the constant is 
just (5ij + j?)- From this Q4.12 ) follows in general. 

Concerning this subalgebra let us make a side remark. Let us define a new set of fields 

Si = S 1} Sj = Sj - Sj-i, 2<j< M. 
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Then the Poisson brackets among these new fields can be expressed in the following matrix form 

/ N+l 



N 
-1 



V o 



-1 


.. 


. 








2 


-1 . . 


. 











. . 


. -1 


2 


-1 





. . 


. 


-1 


2 



(4.14) 



This is almost the Cartan matrix of the sl(M + 1) Lie algebra, except for the element in the upper 
left corner, which is the Poisson bracket {Si, Si}. This may hide some still unknown relation between 
W(N, M) algebras and ordinary finite dimensional Lie algebras. 

(vi). Some useful Poisson brackets 

Proposition 4-5 In the W(N,M) algebras, we find in particular the following brackets: 



{flj,5i} 2 



J_ 

N 



N + l 

i + i 



+ y, 



(N + l)j - NQ, + 1) 
N(N -I) 



N-l 
N-j 



aid j ~ 



+S\ — dependent terms j < N — 1; 

N-l 

W, Si} 2 = ((5 + Sif + Si^-^S'ix - y); 

i=i 

{a N+i , Sj} 2 = 0, i>j; 
{S i ,S j } 2 =(s ij + ±-)6'(x-y). 



(4.15a) 

(4.15b) 

(4.15c) 
(4.15d) 



These brackets are particularly important since they are crucial for Hamiltonian reductions. 

Proof : Let us start from eq.( 4.15c| ): we see that it is true for all the algebras given in Appendix 
B. Now let us consider the case of arbitrary N and M = 2 in which the gauge invariant functions 
have the following realization 

ui = a N (-d + Si) 1 -1+ Y a N+1 (-d + S 2 ) l2 (-d + S 1 ) h -1. 



Comparing with eq. 



we see that 
ui = ui + Sui, 



5ui = on+i — dependent terms 



Since u{s and ui's (together with {a;(l < I < N — 1)}) form the same algebra ( [4.6a -4.6c), we may 
use this invariance to derive the additional Poisson brackets. For example 

{u , u 1 } 2 = {u , ui + aAr + i} 2 

The aN+i -dependent terms on the LHS determine {wO) a JV+i}2 = {&N, ^N+l}2- Similarly 

{ui,ui}2 - {ui,ui} 2 = {aN+l,Ui} 2 + {«i,ajv+i}2 + {«AT+l,OJV+l}2, 

The RHS must be i-dependent, since all the a^+i -independent terms are excluded. This require- 
ment implies that 



{a/v+i) Si} 2 must be a^v+i — dependent 
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However, by dimension counting, it is easy to see that a/v+i -dependent terms are not allowed in local 
polynomial Poisson brackets. Therefore we must have 

{a>N+i,Si}2 = 0. 

For the same reason, when we consider M > 3, i.e. we add more pseudo-differential terms to 
the Lax operator, we will obtain ( 4.15c) ). Therefore, when we add more pseudo-differential operator 



terms to the Lax operator, they will not change the Poisson brackets involving the Si fields already 
considered, thus we must have in general ( 4.15a| ). 



(vii). The generating algebras of the usual W infinity algebras 

Finally let us consider the N = 1 case. As we have pointed out several times the gauge invariant 
functions {u/} form a Wi +oa and a Woo algebra. However these two algebras can be realized also by 
means of 2M fields only, via the combinations Q4.4D and the finite dimensional algebras determined 
by the first and second Poisson brackets of the fields. We could therefore phrase this situation by 
saying that the TVi+oo and the W^ algebras reduce to the latter finite algebras. 

Since M can be any positive integer we see that there are infinite many different realizations of 
the VFi+oo and the Woo algebras. For instance, W(l, 2) and W(l,3) give the four- and six-field 
realizations of the Woo algebra. Alternatively we can say that Wi+oo and Woo algebras are highly 
reducible, and multi-field representations of the KP hierarchy provide a way to classify the reduced 
algebras. 



5 Reductions of the (AT, M)— th KdV hierarchy 

We have claimed above that W(N, M) algebra can be reduced to the usual Wn+m algebra. In this 



section we will show this and the more general relations (|l.5| ) and (1.6). Let us first summarize the 



two schemes for Hamiltonian reduction with second class constraints we will be using. 



5.1 Two reduction schemes for Hamiltonian systems 

Let us consider a Hamiltonian system, no matter whether it is integrable or non-integrable. It will 
have fields generically denoted by fa, a Hamiltonian H and a Poisson bracket. The equations of 
motion are 

^fi = {fi,H}, i = l,2,...,n. 

Now we want to impose, for example, C = where C is a particular combination of the fields and 
their derivatives, and suppose that it is second class, i.e. 

{C,C} = A5(x-y), A| c=0 /0, 

then, in order to avoid inconsistencies, we can proceed in two ways. 

The first scheme 

If we explicitly know all the Poisson brackets among the independent fields, we can first improve 
the Poisson bracket, then derive the reduced equations of motion. 

step 1 introduce in the reduced phase space the Dirac-Poisson bracket 

{fiJjte = ({fiJA - {ft'Cy^icJi}) ■ (5.i) 
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step 2 :The Hamiltonian of the reduced system is just H\f 1= Q. Using this restricted Hamilto- 
nian and the Dirac-Poisson bracket, we are able to derive the equations of motion for the reduced 
Hamiltonian system. 

The second scheme 

If it so happens that we do not know the full Poisson algebra, but we know the Poisson brackets 
between the constraint and all the fields, then we can improve the Hamiltonian with the addition 
of a suitable Lagrange multiplier term, and make use of the known Poisson brackets to derive the 
equations of motion of the reduced system. 

step 1 : Add a Lagrange multiplier term to H 

H => H = H + J aft, 

where a is an expression to be determined. 

step 2 : Using this new Hamiltonian and the original Poisson brackets derive the equation of 
motion for C, 

step 3 : Determine a so that the second class constraint is preserved by the time evolution generated 
by the new Hamiltonian, i.e. turn C = into a first class constraint. 

step 4 '■ Using this new Hamiltonian and the original Poisson brackets derive the equations of 
motion for the reduced system 

d 

Oj.fi = {fi,H}\ c =o, i = 2,3,...,n. 



Comment concerning integrable Hamiltonian systems 

What we said so far is valid for any Hamiltonian system. Now let us turn our attention to the 
reduction of an integrable Hamiltonian system. After implementing either of the above two schemes, 
we have to make sure that the reduced system is still integrable. In other words we should further 
construct its bi-Hamiltonian structure, which is a recursive but lengthy procedure. A shortcut in this 
sense consists in finding its Lax pair representation. 

Let us finally remark that the two schemes presented above lead to the same results, to the extent 
we have been able to produce explicit results for both of them. 



5.2 Reduction of W(N, M) algebras 

We want to impose the constraint S = for some of the S fields of the algebra, and this is second 
class. 

Applying the first reduction scheme to the W(l, 2)-algebra, one can prove the following sequence 
of reductions || 

W(l, 2) W(2, 1) 3^ W(3, 0) = W 3 . 

For the algebras given in Appendix B, we can find another sequence 

W{\, 3) W(2, 2) 3^ W(3, 1) W(4, 0) = W 4 . 
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These two sequences confirm the relation ( |1.6[) . However, since we do not have at hand a compact 
explicit form of the algebra W(N, M) we can verify that with S\ = 0, W(N, M) reduces to W(N + 
1,M — 1) only case by case. To overcome this difficulty we have to resort to the second reduction 
scheme. 



5.3 Reduction of the hierarchy 

Proposition 5.1 : When we impose the constraint S\ = 0, the W(N, M)-algebra reduces to the 
W(N + 1,M - l)-algebra 



W(N, M) 



S!=0 



W(N + 1,M -1) 



(5.2) 



and, simultaneously, the (N, M) hierarchy reduces to the {N + 1, M — 1) hierarchy. 

Proof : As anticipated above, we adopt the second reduction scheme. Our starting point is 
eqs.(3.11a-3.11c), which are generated by the second Hamiltonian H2 = y < Ln > through the 



second Poisson structure. Obviously the constraint Si = is not preserved by the time evolutions, so 
we modify H2 as follows 



H 2 = H 2 + aS u 



Ho 



with a to be determined. This new Hamiltonian generates the following equation of motion for the 
Si field (with respect to the original second Poisson bracket) 



(—Si 

^ dt2 ' improved N 



2 , 



a[ + 2SiS[ - Si + — ^— « + J{Si,a} 2 Si. 



Setting Si = 0, we obtain 



a 



N+l 



ax 



therefore our improved Hamiltonian is 

H 2 = 2 [ (a 2 (x) + (S N - 



1 



iV + 1 



ai(x)Si(x) \dx. 



(5.3) 



Using eq.(4.15a) and (4.15d), we can derive the explicit form of the equations of motion of the reduced 
system 

d n , 2 / N + 1 \ p+i) 

dr 2 ai=ai+2a ^-wTi[ 1+1 ai 



2 44 / N- k 



N + 1 ^ \ l-k 

k=l \ 



(l-k) 

a k a\ 



1<1<N; 



(5.4a) 



dt- 



1+1 , 

-aN+i = a' N+l + 2a' N+l+l + 2a' N+l S i+ i + 2a N+ i S k ^j , 1 < I < M - 1 

k=i 



dt 2 1 N + 



1-1 



-a , 1 + 2^^-S;'-2(^5 fc ) , 2<1<M. 
k=i 



(5.4b) 
(5.4c) 
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We immediately recognize that this is nothing but the first non-trivial flow of the (N + 1, M — l)-th 
KdV hierarchy. The only effect of our reduction is the shift 

N — ► — >N+l, M — ► — >M-1. 

The algebra associated to this hierarchy is W(N + 1, M — 1), thus we are led to the relation ( |5.2| ). 

But now we remark that our proof is valid for any value of M. Therefore we can adapt the above 
proof to the case when we set 52 = 0, £3 = 0, . . . , Sm = 0. Finally we get the (N + M)-th KdV 
hierarchy and the Wn+m algebra. 



6 Drinfeld— Sokholov representation 

The integrable models and reductions we have been considering so far can be synthesized in a very 
compact and useful form via suitable Drinfeld-Sokholov (DS) linear systems. From them one can 
easily extract the corresponding Lax pair. 



6.1 The linear system associated to the (iV, M)— th KdV hierarchy 



The (N, M)-th KdV hierarchy ( |1.3[ ) can be viewed as the consistency condition of the following linear 
system 



L*o(A,t) = Atf (A,t), 

£*o(A,t) = (i>) + tto(A,t) 

where A is the spectral parameter, and ^0 is referred to as Baker-Akhiezer function. In |8) we have 
shown that this linear system naturally appears in multi-matrix models. 

In terms of Baker- Akhiezer function, we can introduce an important ingredient - the r-function, 

w (A,t) = ^ , (6.2) 



where 



OO OO -I r\ 

r=l r=l 



r\ r dt r ' 



is a vertex operator. One can show that the the r-function satisfies the following relations [11] 

o 2 



dtidt 



■lnr(t) = res 9 LN. (6.3) 



r 



6.2 The generalized DS representation 



Now let us return to the spectral problem fl6.1|) . The first equation of ( |6.lD is highly non-linear, our 
aim is to linearize it. In order to do so, we introduce some more notations, 

(-Etf)w = SikSji, l<i,j,k,l<N + M, 
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which is the (N + M) x (N + M) matrix with only one non-zero element at the position Define 



N+M-l 

1+ = E E i,i+1' 
i=l 

* = (*_ M ,^i_ M ,...,^_i,^o,^i,---,^Jv) T , (6.4) 

M N+M-l 

S = E SM-i+lEi,i, A = E a-iEN+M,N+M-i- 

i=l i=l 

So VP = \l/(A,t) is a column vector, the other are matrices, is the Baker- Akhiezer function of 
eq. fl6.ip . Now we can express the spectral equation in linear form as follows 



d + S + A- XE N+M ,M+i -1+)* = 0. (6.5) 

If we eliminate all the VP; in favor of Vl/o, we recover the spectral equation in (|6.1|). 

Now let us consider some particular cases 
(i). M = case, this is the N—th KdV hierarchy, the above linear system is just the original Drinfeld- 
Sokholov representation, in which the spectral parameter A is put on the lower left corner [16|. 



(ii). N = 1 case, this is 2M-field representation of KP hierarchy. The spectral parameter in this 
case is located on the main diagonal line (on the lower right corner). 

For general W and 'M', the spectral parameter appears in the last row. When we make the 
reductions discussed in section 5, Imposing Si = is equivalent to replacing 

N — ►JV+l, M — >M-l. 

In fl6.5|) , this change corresponds to moving spectral parameter A from the right to the left by one 
step. Repeating the procedure, we finally can move A to the lower left corner, that is, we get a KdV 
hierarchy. 



7 The dispersionless (AT, M)— th KdV hierarchy 

In this section we will consider the dispersionless limit of the (N, M)-th KdV hierarchy (|l.3| ). The 
usual procedure to get this limit is the following: we simply ignore the higher than first derivatives 
in the equations of motion as well as in the Poisson brackets. This is equivalent to substituting the 
commutators by the basic canonical Poisson relation, i.e. 

[d,x]=>{p,x} = l, (7.1) 

where p denotes the canonical conjugate momentum. In terms of this elementary Poisson structure, 
the dispersionless (N, M)-th KdV hierarchy can be written as 



gpc = {£;,£}, (7.2) 



with the dispersionless Lax operator 



N M , 

fr[ fri (p - Si)(p - s 2 ) . . . (p - Si) 

N oo 

= p^+E^ + E^" 1 - ( 7 - 3 ) 



i=i i=i 
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The subscript '+' in (7.2) means that we keep only non-negative powers of 'p'. As an example we 
give the second dispersionless flow equations 

^-a; = 2a' l+1 - 2< " N ~ l \ i-ia[, 1 < I < N - 1; (7.4a) 
— a N+l = 2a' N+l+1 + 2a! N+l S l+1 + 2a N+l S k ) , < I < M - 1; (7.4b) 



k=l 



-^St = + 2S,S;, 1 < Z < M. (7.4c) 

Comparing with eqs( |3.11^ - |3.11c ), we see that only the first order derivatives survive. In the disper- 
sionless limit, the gauge invariant functions it;'s have very simple expressions. The first few are 

UQ = CIN, Ul = djV+2 + a>NSl, 

u 2 = a N+3 + a N+2 (Si + S 2 ) + a N Sf, 

In order to find general compact formulas we introduce a set of completely symmetric and homogeneous 
polynomials e k 

ek(Si,S 2 , ... ,Sm) d = ^2 s h s i2 ■ ■ ■ Si k , k G Z + (7.5) 

l<i 1 <i 2 <...<M 

In particular we have 

k 

e k (Si, S 2 , . . . , Sm) = ^2 ^i e k~i(S 2 , S 3 , . . . , Sm)- (7.6) 

1=0 

Without further specification hereafter we understand e& = efc(Si, S 2 , . . . , Sm)- It is not difficult to 
show that these symmetric functions satisfy the following identities (see Appendix D for the proofs) 

d ' 

— e i+ i = S i e i-^ ( 7 - 7a ) 

' /i=0 
M i 

53E' s r e i-Ai = ( i + M ) e i; ( 7 - 7b ) 

1=1 i_i=0 

M q Mi 

e 't+i = Yl ^W e i+i = X) £ BiSfe-v ( 7 - 7c ) 

i=i a ^ 1=1 /x=0 

M i— 1 i j 

E] EI EI S^ei- a e j+a -f3 = ^2(i-a + M)ei- a e j+a ; (7.7d) 

|=l/3=Oa=/3+l a=l 
M i a~l i 

EEE( a "^" l)5|5f ei-cej+a-jg-i = ^ ae-_ Q e i+a ; (7.7e) 

1=1 a=l 13=0 a=l 

El(^7r e i+i) 5 ( a7T e i+i)= (j + M) ei dej +^2(le j+ idei„i + (j - i + l)ei-ide j+ i) . (7.7f) 
l=i oa i ^ i=i 
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The last one is an operatorial equation. Using these symmetric polynomials we obtain 

l+i 

ui = ^2 a N+k-iei-k+i(Si, S 2 , ■ ■ ■ , S k ), < I < M - 1; 
k=i 

M 

u i = ^2 a N+k-iei-k+i(Si,S 2 ,...,S k ), 1>M. (7.8) 
k=i 

In the dispersionless limit the two Hamiltonian structures are extremely simplified. With a little 
exercise, we get the first Poisson algebra 

{cLi,aj}i = (N5 i+j> Nd + (N - z)a i+i _Ar_i<9 + (N - j)da i+j ^ N ^5(x - y), (7.9a) 
{ui, uj}i = {iu i+j ^id + jdu i+j -i)5(x - y), (7.9b) 
{a i ,u j } 1 = 0. (7.9c) 

The second Poisson algebra is 

{ai,aj} 2 = [ia i+j -id + jda i+j -i + iu i+j - N -id + jdu i+j - N -i 
i-l 

+ l)a i+j -i- 2 dai + (j - I - l)ajda i+i _2_ 2 ) (7.10a) 

i=i 

+ ^2((J ~ l ~ l)a>ldiH + j-.i-N-2 + l)u i+j ^ N ^ 2 dai 
i=i 

(N-i){N-j) 
H j-f Oi-idaj-^Six - y), 

{ui,Uj} 2 = [(N + z)n i+i+Ar _i5 + (N + j)5n i+i+A r_i + —u^duj-i 
N-i 

+ Y, (( N + 3 ~ 1 ~ l)oj9w*+i+JV-i-2 + (N + i-l- l)u i+j+N ^ 2 dai) (7.10b) 

i-l 

+ ~ Z ~ l)u i+j+N ^ 2 dui + +(j - / - l)utdu i+j+N -i- 2 ^]5(x - y), 

1=0 

(N-i)j 
N 

i-2 

+ ~ 1 ~ l ) u i+j-i-2 da i + (N + j - I- T)aidu i+j -i- 2 )]S( x ~ v)- 

i=i 

The conserved quantities are [] 

N 



{at, uj} 2 = [iu i+j _id + (N + j)du i+j -i ai-iduj-t (7.10c) 



H r = — J res p (£^y r > 1. (7.11) 



One can check that the two Poisson brackets are compatible w.r.t. these quantities, and they indeed 
generate the flows ( |T.2| ) . In the remaining part of this section we will try to derive the Poisson brackets 
among different Sj's, and the Poisson relations between Si and other fields. 



"Here the residue means simply the coefficient of the p 1 term. 
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Proposition 7.1 The symmetric polynomials satisfy the following Poisson brackets 

'i + M 



{e m ,e i+1 } 2 = \^—— + + M)e t de 



N 

i 

+ ^2{le j+ idei-i + (j - i + l)ei-ide j+ ^ 5(x - y) 
1=1 



(7.12) 



Proof : As we discussed in subsection 4.3 the possible Poisson brackets among the Si fields are 
of the form ( 4.13j ). Taking the Poisson bracket between two monomials of the Si fields reduces the 
total number of Si fields by 2. Eq. ([4.4j) shows that in Ui + u{i > 0) the terms with lowest powers of Si 
fields are a^+M-i e i+l- Now let us consider the Poisson bracket {ui + M, uj + m}2 with i,j > 0. In this 
expression the terms with lowest powers of Si fields will come from 

aN+M-i{ei+i, ej+i}207V+M-i 
On the other hand, the comparable terms in the RHS of the Poisson bracket ( 7.10b| ) are 

[/i + M 



aN+M-i 



N 



+ l) (j + M)e i 8e j + Y,(^j+idei 



i=i 



+(j - i + l)ei-ide j+ i) qn+m-iS(x - y) 



Comparing these two expressions, we obtain (7.12). 

Proposition 7.2 The Poisson brackets ( [7.1 2]) imply the Poisson algebra ( 4.15d ). 

Proof. This is not difficult to prove. We have already noticed that the bracket {Si, Sj} must have 
the form (4.13). Now making use of the relations ( |7.6| ) and (7.7a-7.7f), we can prove that ( |4.15d ) 
indeed leads to eqs.( [< 7 .12" ). This uniquely fixes the undetermined constant of ( 4.13 ). 

The Poisson brackets between S\ and ai{\ < I < N + M — \) are easier to derive. In fact we can 
get them directly from eq.( [4.15a - 4.15d ) by suppressing all higher order derivatives, the result is 



i-2 



K Sih = [S[ + J2 aiSr 1 - 1 + ^a i _ 1 ]«5'(x - y), 
i=i 



{a N , 5x} 2 = (S? + aiS?- l )6'(x - y); 
i=i 

{a N+l „ l ,Si} 2 = 0, l = 2,3,...,M. 



1< i < N - 1: 



(7.13a) 

(7.13b) 
(7.13c) 



8 Reductions of the dispersionless (TV, M)— th KdV hier- 
archy 

Starting from the previous results we will examine in this section the reduction of the dispersionless 
(N, M)-th KdV hierarchy . To find reductions we may suppress the fields Si one by one. This 
was already done in section 5 and in || for the dispersive (general) case and will not be repeated 
here. In this section we are interested in the possible existence of other reductions. Let us take the 
dispersionless four-boson representation of KP hierarchy as an example. The Lax operator is 

p-Si {p-Si){p-S 2 ) 
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The second Poisson algebra can be obtained from (B.l) by killing all the higher derivatives 

{a 1 ,a 1 } 2 = (aid + dai)S(x -y), {a\,a 2 } 2 = (a 2 d + 2da 2 )S(x -y), 
{ai, Si} 2 = Si5'(x - y), {ai, S 2 } 2 = S 2 5'(x - y), 

{a 2 ,a 2 } 2 = [a 2 (2S 2 - S^d + da 2 (2S 2 - Si)]$(x - y), (8.2) 

{a 2 , S 2 } 2 = (at + S 2 (S 2 - Si)^5'(x - y), 

{a 2 ,Si} 2 = 0, {Si,Sj} 2 = (5ij + l)6'(x - y), i,j = 1,2. 

We call it the classical w(l, 2) algebra. With respect to this Poisson algebra, the second Hamiltonian 

H 2 = J (a 2 + a 1 S 1 ), (8.3) 

will generate the first non-trivial flow equations 

d 

— ai = 2(a 2 + aiSi)', 
ot 2 

— a 2 = 2a' 2 S 2 + 2a 2 (Si + S 2 ) , 
ot 2 

jLs 1 = (2a 1 + sfy, 

^S 2 = (2ai + Si)'. 
ot 2 

Now, instead of imposing S\ = as we have done previously, we let S\ = S 2 . Starting from the 
algebra (B^) we get an improved algebra 



{ai,ai} = [aid + da\j5(x -y), {01,02} = (a 2 d + 2da 2 )5(x -y), 

{a 2 , a 2 } = (a 2 Sd + da 2 S - -a\da\)8(x - y), (8.4) 

f 3 
{ai,S} = S5'(x -y), {a 2 ,S} = -axS 1 \x - y), {S, S} = -5'(x - y). 



With respect to this reduced algebra, the original Hamiltonian (3.3) generates the following flow 
equations 

d 

— a\ = 2(a 2 + aiS)', 
ot 2 

^-a 2 = 2a' 2 S + 4a 2 S', 
ot 2 

J-S 2 = (2 ai + S 2 )'. 
ot 2 

It turns out that these equations admit the following Lax representation [] 

d 

^-£ 4 = {(£4)+) A}- 



*This implies integrability for the reduced system. The Si = S2 reduction for the corresponding dispersive hierarchy was 
studied in SJ; the reduced system however is not integrable, at least as long as we stick to locality and polynomiality. 
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with 



C2=p + ^s + Tf^W (8 - 5) 

As we explained in Q, imposing S = will reduce the above hierarchy to the dispersionless Boussinesq 
hierarchy. 

The obvious generalization of the S± = 5*2 reduction to other hierarchies consists in picking out 
i (i = 1,...,M) of the Si fields and setting them equal. So altogether we expect to find k distinct 
reductions with 

-(fHiO--+U)- M -— 

We have checked this for the 6-boson field representation of KP hierarchy. 



9 Conclusions 

The content of our paper can be summarized as follows. We have systematically discussed the (N, M)— 
th KdV hierarchy . The large integrable differential hierarchy (|1.3| ) contains both the higher KdV 
hierarchy and the multi-field representations of KP hierarchy. The second Hamiltonian structure of 
this hierarchy leads to the extended W-algebra, W(N, M). By suppressing the Si fields in succession 
we can reduce the W(N, M) algebra to the usual PFjv+M-algebra. Simultaneously the corresponding 
(iV, M)— th KdV hierarchy reduces to {N + M)-th KdV hierarchy. There do not seem to exist other 
integrable reductions as long as we put restrictions only on the fields Si . However in the dispersionless 
limit there do exist another kind of reduction in which we identify some of the Si fields. The resulting 
integrable hierarchy can be obtained by directly imposing this constraint in the Lax operator. 

In we use the results of this papers to compute correlation functions of the two-matrix models. 
As far as the study of hierarchies and algebras is concerned, there remain some open questions. One is 
whether there are other reductions of the (N, M) hierarchy of different nature than the ones we have 
considered. Another problem concerns the quantum versions of the algebras W(N, M). We have seen 
in section 4 that each algebra W(N, M) contains a Virasoro sub-algebra (|4.7| ); this may suggest the 
quantum version of W(N, M) algebra to correspond to some new conformal field theory. |17| contains 
some hints on this direction: the W{2, 1) algebra is examined and its several free field realizations 
and quantum version are derived. 

A third interesting question is whether our W(N, M) algebra has any relations with the wf^ 



N 



introduced in [18]. The latter is deduced from a WZW model via Hamiltonian reduction. Finally it 



would be quite interesting to understand the relation between the extended KdV hierarchy (1.3) and 
the conformal affine Toda theories (CAT models), for we know the two boson representation of KP 
hierarchynaturally appears in the si (2) CAT model. 



Appendices 
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A Derivation of the second flow equations 

In this Appendix we will derive the second flow equations eqs.( |3.lia| - |3,llc ). First we see that 



dt 2 h^ m ^ £i dt2 ^9- Sid- Si-! "" d-Si 

Our next task is to calculate the commutator [(^Ln^J : L\. Obviously, we have 



After a straightforward calculation we obtain 

Y - J. 

«, . r) N + 

A 



[(it*) , L+ ] = [a2 + ^a 1 ,^+5:a^- 1 ] 



= 2 £ a ^-< + E -r^- 1 - 1 - If £ , <4 r) ^ 

1=1 1=1 r=l \ I 

N-2 N-l-1 / AT , 1 \ 

Z=l r=l V ' / 

2 // 
[d 2 + -^ai, ajv+i-i] = 2a / Ar+ i_ 1 (9 + a^ +i _! 

= 2a^ v+z _ 1 (5 - Si) + 2aj V+ j_ 1 5j + a N+l _ 1} 

[d 2 + loi,a - S k ] = -2S' k d - S' k - loi (A.2) 

= -(a - s^sj; + 5;' - 2s k s' k - 



A d — S k d-S k l A J 5-5 fc 
2S ' k dh- k + (I ' 1 + 2SkS ' k ~ Sk ) d~= 



Sk 



The last two formulas enable us to calculate 



M 



1 1 1 



d-s^-s^ d-Si 



M-l 1 1 1 



2a' N + 2 E a> N+i-i 



d - s t d - s^t ' • • d - S! 
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M '111 
+ E(^-i + 2a^-i^) d _ Sl O-S l ^---d^S- 1 



(A.3) 



M I 



+ § £ aN+l ' l dhi ■ ' ' 9^ (I"' 1 + 2Sfe ^ " 5 d^~ k ■ ■ ■ d^S 1 

M l-l j j 11 

+ E E *N+l-l J^T ■ ■ ■ (24) ^T— 



2=2 fc=l 



d-s t d-s k+1 y K Jd-s k d-Si 



The last term here is not in the standard form: we have to move S' k to the left 
M l-l , j j j 



Z=2 fc=l 



^(2 8w _ 1 S[) — 



1 1 



Z=2 fe=l 
M !-l 



d - s l d - si_! a - 5i 
1 1 „„ 1 



M Z-l 111 
= E(E 2a^_ 1 ^,) a _^ a _^ i . . . 



(A.4) 



Z=2 fc=l 

M i-1 , , fc-1 



jw t — j_ ^ ft; — J. 



Combining eqs,( A,2 ), ( |A.3| ) and (|A.4| ), and comparing with eq(|A.l|), we obtain eqs.( 3.11a - 3.11c| ), 



which are the first non-trivial flow equations in the hierarchy (|1.3| ). 

B Some simple W(N, M) algebras 

In this Appendix we give explicit expressions for a few simple W(N, M) algebras. 
W(l,2) algebra 

{ai,ai} 2 = {aid + dai)5(x -y), {ai,a 2 } 2 = (a 2 d + 2da 2 )5(x -y), 
{ai,5i} 2 = (5 2 + SiSWx - y), {ai, 5 2 } 2 = (2d 2 + S 2 d)5(x - y), 
{a 2 ,a 2 } 2 = [(2a' 2 + 4a 2 5 2 - 2a 2 S 1 )d + 4' + (2a 2 <S 2 - a 2 S x ) r \5(x - y), (B.l) 
{a 2 ,S 2 } 2 = (a l + (d + S 2 )(d + S 2 -S l ))5\x-y), 
{a 2 ,Si} 2 = 0, {Si,Sj} 2 = (5ij + l)S'(x -y), i,j = l, 2. 
The associated Lax operator is 
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This algebra generates the algebra through the transformation (|4.4|) . For this reason it is called 
the 4-boson representation of algebra. 

W(2, 1) algebra 



{ai,ai} = (^d 3 + aid + da^j5(x - y), 



>2 

{o 2 ,a 2 } = (d 2 a 2 -a 2 d 2 + 2a 2 S 2 d + 2da 2 S 2 )8(x -y), (B.3) 
{ai, a 2 } = (a 2 «9 + 2da 2 )<5(* - y), {01, 5 2 } = (h 2 + S 2 d)5(x - y), 
{a 2 ,S 2 } = (ai + (d + S 2 ) 2 )5'(x-y), {S 2 , S 2 } = |<5'(x - y). 
The associated scalar Lax operator is 

L = d 2 + ai + a 2 —^— . (B.4) 
o — b 2 

W(3, 0) = W 3 algebra 

The W(3, 0) algebra is nothing but the W 3 algebra 

{01,01} = (2<9 3 + aid + dai)5(x - y), 

{01,02} = (a2<9 + 2<9a 2 — d 2 ai — d A )S(x — y), (B.5) 



{a 2 , a 2 } = [2a' 2 d + o 2 - ^(01 + d 2 )(dai + d 3 )]<5(x - y). 
The associated Lax operator is 

L = d 3 + aid + a 2 . (B.6) 

W(l,3) algebra 



{01,01} = (aid + dai)5(x - y), {a 1 ,a 2 } = (a 2 d + 2da 2 )S(x - y), 

{01,03} = (a 3 d + 3da 3 )5(x - y), {oi, Si} = (d 2 + S 1 d)5(x - y), 

{ai,S 2 } = (2d 2 + S 2 d)5(x-y), {a u S 3 } = (3d 2 + S 3 d)5(x - y), 

{02,02} = (d 2 a 2 -a 2 d 2 + 2a 3 d + 2da 3 + a 2 (2S 2 -Si)d + da 2 (2S 2 -Si))5(x-y), 

{o 2 , 03} = (3d 2 a 3 - a 3 d 2 - 2Sida 3 - a 3 dSi + 2S 2 da 3 - da 3 S 2 

+ 2a 3 S 3 d + 3da 3 S 3 )5(x-y), (B.7) 

{a 2 ,Si} = 0, {a 2 ,S 2 } = ( ai + (d + S 2 )(d + S 2 -Si))5'(x-y), 

{a 2 ,S 3 } = (ai + 3d 2 + (3S 3 + S 2 -2Si)d + (S 2 -S 1 S 3 -S' 1 -S' 2 + 2S' 3 ))5 , (x-y), 

{a 3 ,a 3 } = (a 3 d 3 + d 3 a 3 + a 3 dai + a x da 3 + a 3 (Si + S 2 - 3S 3 )d 2 
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- 8 2 a 3 (S 1 + S 2 - 3S 3 ) + a 3 (3,S3 2 + 5i5 2 - 2SiS 3 - 2S 2 S 3 + S[ - 3S' 3 )8 
+ 903(351 + S ± S 2 - 2SiS 3 - 25 2 5 3 + S[ - 3S' 3 j) 5'(x - y) 

{a 3 ,S 3 } = (a 1 d + a 2 + a l (S 3 -S 2 ) + (d + S 3 )(d + S 3 -S 1 )(d + S 3 -S 2 ))5'(x-y), 
{a 3l Si} = 0, i = 1,2; Sj} = (<% + 1)6' (x - y), *,j = 1,2,3. 

The Lax operator is 

L = 9 + ai wh; + a2 dh- 2 dh; + ^dhswh^dh;- (R8) 

This algebra also generates the algebra, and is called the 6-field representation of the algebra. 
W(2,2) algebra 



{ai,ai} = (^<9 3 + aid + da^j5(x - y), 

{01,02} = (a 2 d + 2da 2 )5{x - y), {oi,a 3 } = (a 3 8 + 3da 3 )5(x - y), 

{01, S 2 } = {^8 2 + S 2 d)5{x - y), {ai, S 3 } = (^8 2 + S 3 d)6(x - y), 

{a 2 , a 2 } = (d 2 a 2 - a 2 d 2 + 2a 3 d + 2da 3 + 2a 2 S 2 8 + 28a 2 S 2 )s(x - y), 

{a 2 ,S 2 } = {a 1 + (d + S 2 ) 2 )S'(x-y), (B.9) 

{a 2 , S 3 } = ( ai + 39 2 + (S 2 + 3S 3 )8 + (S 2 + 2S' 3 - S' 2 ))5'(x - y), 

{a 2 , a 3 } = (3d 2 a 3 - a 3 8 2 + 2S 2 8a 3 - 8a 3 S 2 + 2a 3 S 3 d + 38a 3 S 3 )5(x - y), 

{a 3 , a 3 } = (a 3 d 3 + 8 3 a 3 + a 3 8 ai + ai 8a 3 + a 3 (S 2 - 3S 3 )8 2 - 8 2 a 3 (S 2 - 3S 3 ) 

+a 3 (3Sl - 2S 2 S 3 - 3S' 3 )8 + 8a 3 (3Sl - 2S 2 S 3 - 3 < S 3 ))o'(x - y) 
{a 3 , S 3 } = {a^ + a 2 + ai (S 3 - S 2 ) + (8 + S 3 ) 2 (8 + S 3 - S 2 ))s'(x - y), 

{a 3 ,5 2 } = 0, {S i ,S j } = (6 ij + ±)6'(x-y), i,j = 2,3. 
In this case the scalar Lax operator is 

L = 92 + 01 + a *dhr 2 + ^whr.dh;- (R10) 

W(3, 1) algebra 



{ax, ax} = (28 3 a 1 8 + dai)5(x - y), 

{oi,a 2 } = (a 2 8 + 2<9a 2 - 8 2 ai - 8 A )5(x - y), 

{01, a 3 } = (a 3 8 + 38a 3 )5(x - y), {a u S 3 } = (28 2 + S 3 8)5(x - y), 

2 

{a 2 , a 2 } = (d 2 a 2 - a 2 8 2 + 2a 3 8 + 2da 3 - -(ai + 3 2 )3(<9 2 + oi))<5(s - y), 
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{a 2 ,a 3 } = (3d 2 a 3 -a 3 d 2 + 2a 3 S 3 d + 3da 3 S 3 )5(x-y), (B.ll) 

{a 2 , S 3 } = (|ai + jj<9 2 + 3S 3 d + (S| + 2 < S 3 ))<5'(x - y), 

{a 3 , a 3 } = (a 3 d 3 + <9 3 a 3 + a 3 dai + ai<9a 3 - SasS^d 2 + 3d 2 a 3 5 3 

+3a 3 (5 3 2 - S' 3 )d + 3da 3 (Sl - S' 3 ))s'(x - y) 
{a 3 , S 3 } = (a 2 + ai(0 + 5 3 ) + (5 + S 3 ) 3 )<5'(x - y), 



{5 3 ,5 3 } = ^ / (x-y). 



with Lax operator 



L = 8 3 + aid + a 2 + a 3 1 - . (B.12) 

a — d 3 



TF(4,0) = W 4 algebra 

At last let us give here the W4 algebra 

{ai, a{\ = (5<9 3 ai<9 + dai)<5(x — y), 

{ai, a 2 } = (a 2 <9 + 2<9a 2 - 2<9 2 ai - 5<9 4 )e)(x - y), 

{ai,a 3 } = (a 3 d + 3da 3 + -{d 5 + 8 3 ai - d 2 a 2 ))d(x - y), 

{a 2 , a 2 } = (d 2 a 2 - a 2 d 2 + 2a 3 d + 2<9a 3 - a^ai - 2ai<9 3 - 2d 3 ai - 6<9 5 )<5(x - y), (B.13) 
{a 2 , a 3 } = (3a 2 a 3 - a 3 9 2 + \{ai+ A8 2 )d{d ai - a 2 + d 3 )),5(x - y), 
{a 3 , a 3 } = (a 3 d 3 + d 3 a 3 + a 3 <9ai + a\da 3 

+-(02 + ai5 + d 3 )d(d 3 + Sai - a 2 )) <5'(x - y). 
It is well-known that this algebra is associated with the scalar Lax operator 

L = d 4 + ai d 2 + a 2 d + o 3 . (B.14) 



C Some VF(iV, M) algebras in the (q, r)-gauge 

W(2, 1) algebra 



{ai,ai} = (i<9 3 + ai<9 + <9ai)<5(x - y), 

{ai, g} = (<?<9 + l;dq)5{x - y), {a 1: r} = (rd + ^dr)S(x - y), 

3 3 
{<1, <l} = ~-qd~ l q8{x - y), {r, r} = --rd~ 1 r5{x - y), 

{q,r} = {f + a 1 + \qdr^r)S{x-y). 
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The associated scalar Lax operator 

L = 8 2 + 01 + qd~ l r. 

W{2,2) algebra 



{ax, ax} = 


(^d 3 + a 1 d + da 1 )5(x-y), 






{ai,n} = 


(n<9+ -drx)S(x - y), {a 1: r 2 } = (r 2 d + 


3 

-dr 2 )S(x 


-y) 


{ai,qi} = 


(qxd + ~dqx)5(x -y), {ax, q%} = (tf2<9 + 


^dq 2 )5(x 


-») 


{ri,rj} = 


~i S ij + ^K<9 _1 r,,-5(x - y), i,j = 1,2, 


{92, rx} 


1 

~ 2' 


{91, n} = 


[d 2 + ai + -qxd' 1 r 1 ^5(x - y), 






{92,^2} = 


(qxfx + —dri + d 2 — drx + l^c^^Wx - 

\ Tx Tx 2 J 


y), 




{91,^2} = 


(—ax + — d 2 + —dr2d + —d 2 r2 — \drxdr 
\rx Tx ri rx rf 


, r2 a 1 

2 H a— 

n ri 




{91,91} = 


/ 3 J 2q 2 „T2 2r 2 „q 2 \ x , s 

-^910 9i + d 1 d—Mx-y) 

\ 2 rx rx rx rx ' 






{91,92} = 


(—d 2 q 2 - —dq 2 d - \qxd' l q 2 - —ax)8(x - 
Vri rx 2 rx ' 


■v), 




{92,92} = 


-^q2d' 1 q 2 )s(x - y). 







In this case the scalar Lax operator is 

L = 8 2 + ax + giS-Vi + q 2 d~ 1 — d'W 

rx 

D Identities satisfied by the polynomials 

We devote this Appendix to proving the identities (0) and (|77a|-^7|. 



Proof of eq-d^ 
le 

G(p;Sx,S 2 . ■ ■ ■ ,S M ) 



Let us define 

I 



(p ~ 5 'i)(p _ ^2) • • • (p - S'm) ' 
which is the generating function of the symmetric polynomials et , 



00 



Taking one time derivative with respect to, say, Sx for G, then expanding it in the powers of p, we 
obtain ( gig ). 
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Likewise we can prove eqs. (|7.7a| - |7.7q j. 
Proof of eq.( 7.7df) 



^ M i-1 i-P 
LHS E E S l e i~P~a e j+& 



1=1 (3=0 a=l 
i i—a 



— E E Sj ei-p-a^j+a 
5=1 (3=0 
i 

= 22 (* ~~ a + M)ei- a ej +a = RHS. 

a=l 

In the last step we have used eq.( 7.7b| ). The proof of eq.( [7.7e| ) is the same as this. 
Proof of eqj7^ 

Since eg. ( |7.7f|) is an operatorial equation, we can cut it into two pieces. 

M 

E 

i=i ~~ i "' 1 /—j 

M 

E 

i=i /—j 

Now let us consider first eq.( D.2a| ) 













( 9 







-j+l 



M i j Mi j+a 

LHS = E E E e-i-aej-pSf +/3 = E E E Ci-aCj+a-pSl 
1=1 a=0 /3=0 l-la=0f3=a 
M i j+a a—1 

= EE(E-E) e <-^+«vf 

1=1 a=0 p=0 f3=0 
i M i-1 i 

= E + a + M )ei-«ej+Q, ~EE E e i-a e j+a-p S l 
a=0 1=1 (3=0 a=(3+l 

i i 

= (j + a + M)ei_ Q ,e i+Q , - ^ (i - a + M)e i _ Q e : , +Q = RHS. 

a=0 a=l 

The LHS of ( p. 2b ) can be divided into two pieces the first being 

M i j-ij-p-i 
A = E E E E S 'kSk S l +Pe i-* e i-P-n-i 

k,l=la=0 (3=0 M=0 
M j—1 M i j+a—1 

= EE^EE E ^r +/3e i-« e i-/3-M-i' 

k=l/j,=0 1=1 a=0 p= a 



Splitting one of the above summations as follows 



j+a— 1 j+a—1 a—1 

E = E 

P=a p=0 p=0 

38 



and using eqs,( 7.7c , 7.7d ), we get 

M j—1 i i 

A = ^^(^(j + a-^-l + M)-^(i-Q + M))e 4 _ Q e J+ „_ At _i5^^ 



a=l 



i M j-l 

= U + M)e i P' j + -i + a)ei- a e j+a - ^ ^ (/j, + l)e i e :/ _ At _i^^ 

a=l fe=l^=0 
M i j-l j+a— 1 

fc=la=l /3=0 



The other term in (D.2b) is 
B 



M i j 

EEE ^i- a ej-pSl + ^ x S' k 

k=l a=0 0=0 

M j-l M i j+a-1 

^^(/J + l)eie J+1 SfS[ + ^^ J2 (0- a + 1 ) e i-a e j+a-/3-iS^S' k . 

k=l/3=0 k=la=l /3=a 



Combining ( p.3[ ) and (D.4) and making use of eq. flT^ej ), we obtain 

i 

A + B = (J + M)e i e' j + ~ i + oi)ei- a e' j+a 



a=l 



M i a-1 

+ EEE( a "^" l)e i _ a e j+a _ /3 _i5f S' k 

k=la=l/3=0 
M i j+a-1 



~ EE E U - 1 + 2a - P - l)ei-aej+ a -(3-lS%S' 
k=l a=l p=j 

i 

= (j + M)eie'j + {ua'i-aZj+a + (j -i + a)ei- a e' j+a 

a=l 

The last term in the intermediate expression vanishes as one can see in the following way. Let us split 
the term into two parts 



(D.3) 



(D.4) 



M i j+a-1 



M i j+a-1 



EE E + a - P - ^e-i-aZj+a-p-iSl S' k , EE E (a-i)ei- a e j+ a-(3-iS^S' k , 

k=la=l J3=j k=la=l /3=j 

and redefine the summation parameter for the first part as follows 

a = i — j — a + /3 + l, f3 = fl, 

After changing the order of the summations the first part it is exactly the same as the second except 
for the different sign. This completes our proof of eq.( p.2b" ). 
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